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Abstract. The main objective of this paper is to describe the dynamic transition of the 
incompressible MHD equations in a three dimensional (3D) rectangular domain from a per- 
spective of pattern formation. First it is shown that the system always undergoes a dynamical 
transition as the Rayleigh number crosses the first critical threshold which corresponds to the 
first eigenvalues of the linearized problem around the basic state. The type of transitions is 
determined by the nonlinear interactions. 

When the magnetic Prandtl number p2 > 1 or the Chandrasekhar number Q < Qq where 
Qo depends on the system parameters, the first eigenvalues are always real. Generically, the 
first eigenvalue is simple and a Type-I or a Type-II transition is possible depending on a non- 
dimensional number exactly given in terms of the system parameters. In particular, the rolls 
are stable patterns after the transition if the Chandrasekhar number Q > or the magnetic 
Prandtl number p2 > p» where Q* and p» depends on the length scales of the domain. When 
the horizontal length scales Li and L2 satisfy •/SkLi = jL2 for some positive integers j 
and fc, two modes, characterizing a hexagonal pattern, become unstable simultaneously. In 
this case, we show that all three types of transition are possible. Depending on two non- 
dimensional parameters, we show that there are eight different transition scenarios. In this 
case we also show that for p2 > 8, the transition is always Type-I with rolls and rectangles 
as stable patterns and hexagons as unstable patterns after the transition. However P2 > 8 is 
a crude estimate and our numerical investigation suggests that this type of Type-I transition 
will be preferred for P2 > P2 where pj < 2.24. 

In the case where p2 < 1 and Q > Qo, the first eigenvalues are complex, and the system 
undergoes a dynamical transition to spatiotemporal oscillatory states. In particular, when 
the first critical modes have time periodic roll structure, we show that the system undergoes a 
Type-I or Type-II transition. In the large Chandrasekhar number limit or the small oscillation 
frequency limit, the transition is always Type-I and the time periodic rolls are stable after 
the transition. 



1. Introduction 

The Rayleigh-Benard convection is a fundamental problem of natural convective heat transfer 
which is characterized by a vertical temperature gradient aligned with the acceleration of gravity 
being maintained over a horizontal layer of fluid. Due to the thermal expansion, the fluid is 
heavier at the top and lighter at the bottom. As the temperature difference between the lower 
and upper fluid boundary exceeds a critical level, a convective motion sets in. There have been 
numerous studies concerning the stability and instability of the convective flows. A detailed 
discussion can be found in [31 131 [3 [B] . 

External magnetic fields change the characteristics of this convection significantly for electri- 
cally well conducting fluids. First it is well known that the critical Rayleigh number and the 
wave number increase with an increasing Chandrasekhar number Q for the onset of convection. 

Date: January 20, 2013. 

1991 Mathematics Subject Classification. 76W05, 35Q35, 35B36. 

Key words and phrases. Magnetohydrodynamic Convection, dynamic transition, pattern formation, hexago- 
nal pattern, metastability. 

The work was supported in part by the Office of Naval Research and by the National Science Foundation. 



2 



SENGUL AND WANG 



Physically this is due to the fact that the energy released by the buoyancy force acting on the 
fluid must balance the energy dissipated by not only the viscosity but also the Joule heating. 
Thus the magnetic field imposes stability to the fluid. Second, the existence of a magnetic field 
allows both the steady and oscillatory convections; see [2 [3]. 

We aim to describe the dynamic stability and transition of the magnetic convection for an 
incompressible fluid in a rectangular domain in M.^ from a pattern formation perspective. As 
is well known, for the magnetohydrodynamics (MHD) equations, due to non-selfadjoint linear 
operator, the transition can be caused by a finite set of real or complex eigenvalues crossing 
the imaginary axis. This makes transitions to both spatially periodic time independent states 
and to spatiotemporally periodic states possible. We focus on the formation of patterns having 
roll, rectangular and hexagonal structures. Our main goal is to precisely determine the type of 
transitions associated with these patterns and hence characterize the stability of these patterns 
in terms of the parameters of the system. 

When the first eigenvalue is real and simple, the transition can only be Type-I or Type-II 
depending on a number exactly given in terms of the system parameters. In particular, when the 
first critical eigenmode has a roll structure, the type of transition is independent of the Prandtl 
number pi. The transition, in this case, is always Type-I if the Chandrasekhar number Q > 
or the magnetic Prandtl number p2 > p* where and p* depends on the length scales of the 
domain. We find that < 307 and p* < 2.24 regardless of the length scales of the domain and 
— > An^ and p* — > 2/\/3 as max{Li, L2} -> 00 where Li and L2 are the horizontal length 
scales. 

Next we study the case where there are two critical real eigenvalues. In this case we only 
consider the special geometry 

Li ^ J 

with positive integers j, k and Li, L2 denoting the horizontal length scales of the box. With 
this assumption, it is possible that two modes which can characterize a hexagon pattern become 
unstable at the same critical parameter. In this case we find that all types of transitions are 
possible in a total of eight different transition scenarios. However, in our numerical investigation, 
we encountered only two of these scenarios. Fixing Q < Q^, we found that the system moves 
from a Type-Ill transition regime to a Type-I regime as p2 crosses p*. In the Type-I transition 
regime, the transition state is an attractor homeomorphic to a circle containing eight steady 
states and the connecting heteroclinic orbits. Of these eight steady states, those having roll 
and rectangular patterns are stable whereas the ones with hexagonal pattern are saddles. In 
the Type-Ill transition regime, a neighborhood of the basic solution is split into two sectorial 
regions Ui and U2- In C/i, the initial conditions move out of this neighborhood whereas in 
i72, the initial conditions will move to an attractor which contains a rectangular pattern as a 
minimal attractor. In this case, we can single out the transition scenario if p2 > 8. We prove 
that when p2 > 8, the transition is always Type-I with rolls and rectangles as stable states while 
hexagons as unstable states after the transition. However p2 > 8 is a crude estimate and our 
numerical investigation suggests that this type of Type-I transition will be preferred for p2 > P2 
where pj < 2.24. It is worth mentioning that when the length scales are in close vicinity of the 
above relation, generically only one of the eigenvalues that characterize a hexagonal structure 
will be unstable. However a second negative eigenvalue will be very close to zero at the critical 
parameter and the structure of transition described above case will still persist. 

Finally, we consider the case where a pair of complex eigenvalues become unstable simulta- 
neously. In this case, we only consider a roll type critical eigenmode. We show that the first 
transition can be Type-I or Type-II. In particular for Q sufficiently large or for the oscillation 
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frequency p sufficiently small the transition is Type-1 and the transition structure is a time 
periodic roll pattern. 

There have been extensive studies on the MHD convection problem, see among others [H [5] 
and the references therein. There are several features of this work that distinguishes it from the 
other studies. 

First, to our knowledge none of the previous studies deals with nonlinear 3D equations in 
a finite box. As we will discuss, the presence of side walls brings in additional difficulties and 
interesting dynamical features. From a physical point of view the addition of side walls is crucial 
especially when the horizontal scales are comparable to the depth of the layer. 

Second, the problem is studied from a transition point of view instead of bifurcation point of 
view. The key philosophy of the dynamic transition theory which is recently developed by Ma 
& Wang [SI [7] is to search for the full set of transition states, giving a complete characterization 
on stability and transition. The set of transition states is represented by a local attractor 
rather than some steady states or periodic solutions or other types of orbits as part of this local 
attractor. 

Third, a crucial step in describing the transition of the system is the reduction of the problem 
on the center manifold. In this paper we use a slightly different method for the reduction of 
the problem which is worth mentioning here. Usually, one expands the center manifold in terms 
of the eigenf unctions of the linear operator of the system. This allows the utilization of the 
approximations for the center manifold in a natural way. However, as in the case of the mag- 
netohydrodynamics convection, this approximation can be difficult to explicitly calculate when 
the eigenfunctions of the linear operator of the original system is complicated for computational 
purposes. We overcome this difficulty by defining a new set of basis for the functional space of 
the problem which behaves nicely under the original linear operator. 

The paper is organized as follows, in Section 2, the mathematical setting of the problem is 
introduced. In Section 3, linear theory is summarized. In Section 4 we state the main theorems. 
We give the proofs of the main theorems in Section 5. In Section 6, physical remarks and 
conclusions are discussed. 



2. MATHEMATICAL SETTING 

We consider thermally driven convection of an electrically conducting fluid in the presence of 
a magnetic field in a rectangular domain = (0, li) x (0, ^2) x (0, h) in M"^. Subject to Boussinesq 
approximation (see among others), the evolution equations read: 

Q%i 1 

— + (u • V)w = (Vp + pgk) + poiW X H) X H + uAu, 

at Po 

fill 

(2.1) -g^ + iu-V)H^{H-V)u + 7jAH, 

dT 

— + (u • V) T = kAT, 
divM = diviJ ~ 0. 

Here u = (ui,U2,u^) is the velocity field, H is the magnetic field, T is the temperature, k — 
(0,0,1), g is the gravitational acceleration, v is the kinematic viscosity, fiQ is the magnetic 
permeability, rj is the magnetic diffusivity (also called resistivity), k is the thermal diffusivity. 
To is the reference temperature at — 0, pa is the density at Tq, and a > is the coefficient of 
thermal expansion. The fiuid density p is given by the equation of state: 



(2.2) 



p = po[l - a{T - To)]. 
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Although the case where the imposed magnetic field and the gravitational field act in different 
directions is also interesting, for simplicity, we will assume that they are parallel. Denoting the 
temperature at = h by Ti, the basic state of (2.1 ) is a motionless state given by: 

u = 0, 

h 

P = Pn- Po9 + a (To - Ti) || 

To put the equations into non-dimensional form, we consider the deviation of the solution 
from the basic state: 

u" u, T" = T - T, H" = H - H, p" = p - p, 

and we set 



u — —u 



H" = '^HnH', T" = (To - Ti) T', 



h rj 



K 



Also we define the following nondimensional numbers: 



pi = v/k the Prandtl number, 

p2 = tj/k the magnetic Prandtl number. 



(2.3) ^^.9a(ro-ri)^3 



the Rayleigh number. 



Using the identity 



fiQ — - — the Chandrasekhar number. 



(V X ff) X = -iy \Hf + (i? • V) H, 



omitting the primes and denoting all the terms that can be written as gradients by p, the 
equations (2.1|-(2.2) take the form: 

'^ + {u-V)u^ pi (-\/p + RT~t + Au + + ^ (F • V) JjV 

ot \ dx3 p2 J 

(2.4) f + + 

dT 

divu = divi? = 0. 



The non-dimensional domain is = (0, Li) x (0,^2) x (0,1) where Li = li/h, L2 = h/h. 
We use the following idealized boundary conditions which are the free-slip boundary conditions 
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for the velocity together with the condition that H remains vertical at X3 = 0, 1. 

dT dH2 dHs 

- — = Hi = — — = — — =0 at xi = 0, Li, 



(2.5) 





du2 _ 




Ui = 


dxi 


dxi 


dui 






dx2 


= U2 = 


dx2 


dui 


du2 




dx3 


dx3 


= "3 



uxi dxi dxi 

1^ — = ^ — = H2 = ^ — =0 at a;2 = 0, L2, 
0x2 0x2 0x2 

T = Hi = H2 = ^=0 at X3 = 0, 1. 
0x3 



We recall here the functional setting of (2.4)-(2.5) and refer the interested readers to [TO 



H = {^iP^ {u, H, T) e [nf I divit = divi? = u-n |ao= H ■ n o} 
i?i = = {u,H,T) e nH \ij satisfies . 

Let Lr = A + Br : Hi ^ H and G : Hi ^ H he defined by 



Ai, = r [piAu + piQ^,P2-^+p2AH,AT] , 
0x3 0x3 



(2.6) 



BRij^V[piRTk,0,u3'j , 
G{'il))^V {H-\/)H~{u- V)u, (i? • V) M - {u • V) iJ, - (u • V) T 



Here V : {Vt) H is the Leray projector. For simplicity, we also use G to represent the 
following trilinear form for -ipi — {ui,Ti, Hi) G H, 



(2.7) 



G(Vl,^2,^3)= / \^{Hi-W)H2-{ui-W)u2 
(i/l • V)U2 - {ui ■ V) H2 



■ U3 



H3 - {ui ■ V) T2 T3 



We also define 
(2.8) 



Gs(V'l,'02,'03) = Gilpl,1p2,ij3) + G(V'2,'01,'/'3)- 

Now (2.4)-(2.5l can be written in the following functional form: 
(2.9) H 



LriJj + G (V) 



3. Principle of Exchange of Stabilities 

The linear stability is determined by the critical crossing of the first eigenvalues from the 
imaginary axis which is called the principle of exchange of stabilities (PES). The linear theory 



associated with (2.4)-(2.5) is well known and can be found for example in [2]. Here we will 
briefiy recall the necessary results. 

First we consider the following eigenvalue problem: 

(3.1) LriP = I3{R)^, 

where Lr is as defined in ( |2.6[ ). Let us denote the set of admissible indices by 
Z = {(Ji,i2, ja) eZ^:j,k> 0, jf+ji ^ 0, ja > l} ■ 
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It is well-known that the eigenvalues of (3.1) having indices J 
following polynomial equation: 



0'i,i2,j3) e Z satisfy the 



(3.2) 
where 



0" ^hifi" ^-hi{R){i + h■l{R)^Q, 



(pi + p2 + 1) 7,7 

P2 



bi =Pi[(p2 + l 



pi 



)7j + P2(3 (J37r)' 



7j 



bi 



Here 
(3.3) 



Pip2(7j + Q(i3^)'7}-a>^, 



a}-(j?ir'+j2%-V 

2 2 I -2 2 

7j = "J +.?3 7r • 



As the linear operator (3.1) is not self-adjoint, the onset of the instability may correspond to 



either real or complex eigenvalues /3, leading to transitions to states with cither time independent 
structures or time periodic structures. We will denote the critical Raylcigh number by Rj. when 



the first eigenvalues are real, and by Rc when the first eigenvalues are complex. Solving (3.2 1 
for R when P — Q and i?e/3 = and minimizing over all admissible indices, we can find the 
expressions for Rr and Rc as follows: 



(3.4) 
(3.5) 



Rr = min 



7j 



p . (P2 + I)(pl+P2)7}r 4 

Rc — mm 

JdZ pi ay •' 



P2pi 



(P2 + 1) (Pl + 1) 



'(J37r)' 



The difference between an unbounded domain and a bounded one is that in the unbounded 



case the minimums in (3.4) and (3.5) are taken over all aj > whereas in the bounded case aj 



has the specific form given above. Figure [T] shows the minimizing indices of Rr as a function of 
the length scales L\ and L^. 



We note that for each J & Z there are three eigenvalues of (3.2) which we can order as 



/3j > /3j > The next proposition states that the principle of exchange of stabilities (PES) 
condition is valid. 

Proposition 1. [2, For p2 < 1, there exists < Qq < oo depending on pi and p2 such that 

2P2(P1 +1) 



(3.6) 



Pl(l"p2)' 



The following assertions hold true: 

//p2 > 1 or Q < Qo then Rr < Rc and Rr is the first critical Rayleigh number. That is. 



there is a finite set of critical indices C minimizing (3.4) such that 

( >OifR<Rr, 

p]{R)\ =otfR = Rr, yjec, 

^-^■'^ [ <QlfR>Rr 

RePj {Rr) < VJ ^ C. 
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&J // p2 < 1 and Q > Qq then Rc < Rr and Rc is the first critical Rayleigh number. That is, 



there is a finite set of critical indices C minimizing (3.5) such that 
(3.8) 



>OifR<Rc, 
Rel3){R) = Rel3j{R) { ^ if R = R^, eC, 

<OifR>Rc 
RePj (Rc) < VJ ^ C. 



There is no simple formula for the critical Chandrasekhar number Qq which separates the 
stationary and oscillatory convection regimes. Qq is determined by either the vanishing of the 
oscillation frequency p (the imaginary part of the first critical complex eigenvalues at i? = Rc) 
given by 

pi Pi + t 

Utilizing > 0, one establishes the estimate (3.6 1. Also it is worth mentioning that in [T], the 
authors prove that Q0P2 > t^^- 

4. Dynamic Transitions 

With the PES condition at our disposal, the dynamic transition theorem in [7] ensures that 
the system always undergoes a dynamic transition to one of the three types (Type-I, Type-II or 
Type-Ill) as the Rayleigh number crosses the first critical value Rr or Rc. The type of transition 
and associated pattern formation are then determined by studying the nonlinear interactions of 
the problem. 



We address in this section the dynamic transition and the pattern formation of (2.4)-(2.5) 



Our main results will depend on two nondimensional numbers a and b which are defined as 
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follows: 

where J = (ji, j2, 1) denotes the critical index and k is given by: 
(4.2) 



7si,S2,2 7si,S2,2 



V P2 ' pi J ' Pl(i?.i,S2,2 - i?r)a} 

2 „,2 



Rs,,s2.s, = %^(7l.„.3 + Q(53^)'), Rr = mini?,„,„,3 = ^(7^ + Qvr^). 



a2 

Si,S2 



4.1. Transitions from simple real eigenvalues. In the case where the magnetic Prandtl 
number p2 > 1 or the Chandrasekhar number Q < Qo where Qq depends on the system 
parameters, the first eigenvalues are real. Generically speaking the first eigenvalue is simple 
and the dynamical transition is characterized by the following theorem. 



Theorem 1. // p2 > 1 or Q < Qq then the problem {2.4) with (2.5) undergoes a dynamic 
transition at R ^ Rr. Consider the case where the first critical eigenvalue is simple. 

a) When the first critical wave index is of the form J = (0,j, 1) or J = (j, 0, 1) with 
non-zero j , the first transition is, 

Type-I, ifb<0, 
Type-IL if b > 0. 

b) When the first critical wave index is of the form J = (j, fc, 1) with j ^ 0, fc 7^ 0. Then 
the first transition is, 

Type-I, if a < 0, 
Type- II, if a > 0. 

c) Moreover, when the transition is Type-I, then as R crosses Rr, the system bifurcates to 
two steady state solutions 'if>± which can be expressed as 

V'±(i?) = ±C{R - Rrf'^i:.] + o((i? - Rrf'^), 

where c > is a constant, and ipj is the critical eigenvector. Furthermore, there is an 
open setU e H with V' = G C/ such thatU = V+UUZ, U+CiU- = 0, V = £ dU+CidU- 
and 'ip± (R) attracts U± . 

The following corollary of Theorem [l] gives sufficient conditions for the existence of Type-I 
transition when the critical mode is a roll pattern. 

Corollary 1. Under the assumptions of Theorem^ when the first critical index is of the form 
J ~ (0, j, 1) or J = (j, 0, 1) with non-zero j , the transition at R = Rr is always Type-I if Q > Q^, 
or p2 > p* where and p* depends on the length scales of the domain. 

a) Regardless of the length scales Li, L2 of the domain, < 307 and p* < 2.24. 

b) As max{Li, L2} 00, we have — )• 47r2 and p* — )■ 2/\/3. 
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4.2. Transitions from real eigenvalues with multiplicity two. In this section we consider 
the transitions when two modes characterizing a hexagon becomes unstable at the same critical 
parameter. For this to happen, a certain relation has to hold for the length scales of the box: 

where j and k are two positive integers; see also Q for a discussion. With this relation in 
hand, the modes with indices / — (j, k, 1) and J — (0, 2k, 1) have the same wave number, 
aj = aj. Moreover, the vector field xjipj + xjipj has a hexagonal pattern parallel to y-axis 
when xi — ±2a;j; see Figure [2j Thus a pair of modes will be critical at the critical Rayleigh 



number whenever one of the modes minimizes the relation (3.4). In Figure 12 the red lines 



show when this is the case at Q = 10. We remark here that when the length scales are in 



close vicinity of the relation (4.3), although the first critical eigenvalue may be simple, the two 
eigenvalues will be very close to each other at the critical parameter Rr and the structure of 
transition will still persist. 

Theorem 2. Consider the case where I = (j, k, 1) and J = (0, 2fc, 1) are the first critical indices 



which satisfy the relation (4.3). Consider the eight regions in the a-b plane as shown in given by 
Figure^with a and b as given in (4.1). //p2 ^ 1 or Q < Qq then the first dynamic transition 



occurs at R = Rr- Moreover the following statements hold true. 

i) When (a, 6) is in one of the regions Ti, T2, the transition is Type-I. The topological 
structure of the transition is given in Figure^and Figure^respectively. In both regions, 
there is an attractor bifurcating on R > Rr which is homeomorphic to containing eight 
bifurcated points and the heteroclinic orbits containing them. Four of these points are 
attractors and four of them are saddles. In region Ti , hexagons 2(f>i + (f>j are the minimal 
attractors, and in region T2, rolls <j)j and rectangles (pi are the minimal attractors. 

ii) When {a,b) is in one of the regions TTi, TZ2, H3, T-Xa^ andTT^, the transition is Type- 
II. The topological structure of the transition is given in Figure \^Figure \ 1 0| respectively. 
In regions TT2 andXT^, there is a repeller bifurcating on R < Rr which is homeomorphic 
to and no bifurcated states on R > Rr. In the regions Hi, H4, TT-z there are 
bifurcated points on both R > Rr and R < Rr. 



iii) When {a,b) is in region XXX\ the transition is Type-Ill; see Figure 11 There is a 



neighborhood 14 of tp = in H such that for any Rr < R < Rr + e with some e > 0, U 
can be decomposed into two open sets Ui , , 



such that 



lim limsupllS'flft, = eU{ 



(4.4) 



limsup||S'fl(t,?/')||H > (5 > VV'GZ^a , 

f 00 

for some 5 > Q. Here Sr is the evolution of the solution with initial data ijj. Moreover 
'P{Ui ') is a small perturbation of the union of two sectorial regions which are given by: 

V{U^^+) KiUf^{xeM?\ -tan^i 1/2 < arg{x) < tan^^ 1/2}, 

V{U^^_) KiUC^{x eW^ \t: - tan^i 1/2 < arg{x) < tt + tan^^ 1/2}, 

where V is the projection onto (pi, (pj plane. 

Moreover, the system bifurcates to two attractors with basin of attraction Uf^^, 
Each attractor consists of three steady states and the connecting heteroclinic orbits. The 
steady states 4>i having rectangular patterns are the minimal attractors of the bifurcated 
attractor. 
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0.0 0,5 1.0 1.5 2.0 2.5 3.0 3.5 



Figure 2. The flow structure of 202.1 + 0o,2 at z = 1 for the box dunensions 
L\ — L2 — 3.5 




a=0 



7^ 



nil 

b=0 


a=0 / 

/ ' 
113 

-- 


-yf 

' /lis 

/ 


114 



(a) 



(b) 



Figure 3. a) Regions of different type of transitions in the coefficient a-b plane 
and b) The topological structure of the solutions in different regions are shown 
in Figure [6] through Figure 1 1 



Corollary 2. Under the assumptions of Theorem^ when p2 > 8, the transition at R — is 
always Type-I. Moreover, there is an attractor bifurcating on R > Rr, which is homeomorphic 
to . This attractor contains eight steady states and the orbits connecting them; see Figure 
^ Two of the steady states have roll patterns, two of them have rectangular patterns and four 
of them have hexagonal patterns. The rolls and rectangles are the minimal attractors of this 
attractor while the hexagons are unstable states after the transition. 



Remark 1. To illustrate the conclusions of Theorem^ we will consider two cases. Solving b 
for p2, we find that b <0 (resp. b > 0) if p\ > cfroii (resp. p\ < Groii) where 

(4.5) Groll ^ 
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(a) R< Rr (h) R> Rr 



Figure 4. Region Ii. 




(a) R< Rr (h) R> Rr 



Figures. Region l2- 




(a) R< Rr (b) -R = -Rr (c) R> Rr 



Figure 6. Region Hi. 




(a) R< Rr (b) -R > Rr 



Figure 7. Region Zl2- 



SENGUL AND WANG 




(a) R< Rr (b) -R > Rr 

Figure 8. Region H3. 




(a) R< Rr (b) -R = -Rr (c) R> Rr 



Figure 9. Region U4. 




(a) R< Rr (b) -R = -Rr (c) R> Rr 

Figure 10. Region H5. 




(a) R< Rr (b) R. = Rr (c) R > Rr 

Figure 1 1 . Region llli . 
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^ 2.0 




Figure 12. The bold lines indicate where two modes satisfying (4.3 ) are critical 
at Q = 10. 



First we will consider the case of small length scales Li, L2 such that the conditions (4.3 I is 



satisfied. Figure 12 shows some of the length scales for which this happens at Q = 10. 
As an example we choose 

(4.6) Q = 10, ii-3/2, L2-\/3Li, 

Under these conditions, two modes satisfying ( |4.3| will be critical with indices / = (1, 1, 1) and 
J = (0,2,1). The critical wave number is aj — 271 / L2 ~ 2.42. In this case y/aroU ~ 0.5. 
Calculating a/h, we find that under the conditions (4.6), the system is in region TZX\ for p2 < 
yJOroU cind will move into I2 o,s p2 crosses yjGroii ■ 

Next we consider the limit of large length scales Li » 1, L2 = ^/iLl. In this case, the 
critical wave number is given by (5.281. If we take Q = 10 again, then we find aj = 2.59. Once 
again, the system moves from region XXX\ to region X2, this time at p2 = y/OroU ~ 0.39. 

4.3. Transitions from complex simple eigenvalues. Under the conditions p2 < 1 and 
Q > Qq, the first critical eigenvalues will be a pair of complex numbers. To demonstrate the 
main ideas and for simplicity, we only consider the transition from a simple complex eigenvalue 
with index Jc — {jc, 0, 1) or Jc = (0, kc, 1) for jc, k^ > 1. 



Theorem 3. Consider b which is defined by (5.53). Assume p2 < 1 and Q > Qq and assume 



that the critical index Jc — {jc, 0, 1) or Jc ~ (0, kc, 1) then we have the following assertions: 
(1) If b < then the problem undergoes a Type-I transition at R = Rc and bifurcates to a 
periodic solution on R > Rc which is an attractor, the periodic solution can be expressed 
as 



4> 



x{t) = 



vit) 



x (t) Re^j^ 
"(A) 



y{t) Imijj^+o(\a{X)\ 



1/2 



a (A) 



1/2 



cos pt, 



1/2 



sin pt. 
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In particular, if Q is sufficiently large or p is sufficiently small then the transition is 
Type-I. 

(2) If b > then the transition is a jump transition and the system bifurcates on R < Rc to 
a unique unstable periodic orbit. 



5. Proof of the Main Theorems 

5.1. Reduction strategy. In dynamic transition problems, the method of proof relies heavily 
on the reduction of the problem to the center manifold in the first unstable eigendirections. The 
center manifold is usually expanded by the eigenniodes of the linear operator of the problem 
and the coefficients of the the expansion are found by the approximation formulas derived in [7] . 
However, for our problem, such an analysis is not suitable because some eigenfunctions of the 
linear operator that are needed to compute the nonlinear interactions have complicated forms. 
In particular the eigenvalues of (3.1) corresponding to indices {j,k,l) with Z > 2, + fc^ 7^ 



can only be found by first computing the roots of (3.2). It is easy to see that these roots will 
have complicated forms involving the parameters of the system. This makes it very difficult to 
compute their contributions in the nonlinear interactions. We will overcome this difficulty by 
expanding the center manifold in terms of the eigenfunctions of the Laplacian operator with the 
same boundary conditions. 

In the case where the first critical eigenvalues are real, the critical eigenfunctions can be 



found by setting (3 = in (3.11. Thanks to the boundary conditions, the eigenvectors can be 
expressed using the separation of variables. For example the vertical component of the velocity 
field becomes: 

(5.1) w=Wjcos — - — cos — - — sinj37rx3, 

Li L2 

and analogous expressions hold for the other components of the solutions of the linear equations. 
Here J = (ji, j2,i3) with nonnegative integers ii,j2, js- It is clear that critical indices must 



have j3 = 1 as the minimum in (3.4) is achieved for J3 = 1 and moreover we must have aj > 0. 
When aj ^ 0, the amplitudes Uj, Vj of the horizontal velocity field and Hi j, i?2.j of the 
horizontal magnetic field depend on the vertical components Wj and H^ j as follows: 



UJ = - - - - 

(5.2) 



Uj = -- — jhTrWj, Vj = -- — ^jsnWj, 
Littj L2aj 



rr ■ TT rr J27r . 
Hi, J = 2"J37rii3.J, -H2,J = 7 2".?3'^-"3,J- 

I^ictj L2aj 

Thus the critical eigenfunctions are determined by three quantities: the amplitude Wj of the 
vertical velocity field, the amplitude H^ j of the vertical magnetic field and the amplitude Oj 
of the temperature field. These can be found as: 

Wj = 7,7, ej = 1, i/3,J - ^- 

Next we need to analyze the conjugate critical eigenfunctions which are the solutions of the 
conjugate eigenvalue problem at /3 = 0: 

Pi{Au* -Vp*)-P2^+T* k =0, 
0x3 

AT* + RpiuX = 0, 
du* 

P2AH* - Qpi— = 0, 
0x3 

divu* = divH* = 0. 
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The conjugate eigenfunctions can be found as: 

W*j = P27j, e} = P2pli?r - -PigTT. 

Now we turn to eigenvalue problem of the Laplacian operator: 

Aej = -7jej. 



The Laplacian operator with boundary conditions (2.51 has a complete set of orthogonal eigen- 
vectors ej — {uj, Tj, Hj) which can be expressed by the same separation of variables ( 5.1 ): 

Again there are only three independent unknowns, namely the amplitudes of the vertical 
velocity w, the temperature T and the vertical magnetic field which are denoted by W , Q 
and respectively. 

• If il + j1 = and j3 ^ 0, then 

(5.4) ej = (0,rj,0), 6,7 = 1. 

• If il +32 7^0 and = 0, then 

(5.5) ej = (0,0, i?j), H3.J = 1. 

In this case, although there are eigenmodes which have the form e = (w, 0, 0) with 
u = {u,v,0), it can be checked that these modes do not affect the calculations of the 
nonlinear interactions as they will be not be present in the lowest order approximation 
of the center manifold function. 

• If jf + j'l 7^ and j3 ^ 0, then the multiplicity of an eigenvalue is three and the 
eigenvectors are: 

e'j = {uj,0,0), Wj = l, 

(5.6) e2=(0,rj,0), 6,7 = 1, 

ei^ = (0,0, i?j), H3,j = l. 

Proposition [T] gives a critical set of indices C = {Ji, J2, ■ ■ ■ , Jk} such that the eigenvalues 
having these indices are the first critical eigenvalues. The dynamics on the center manifold can 
be captured as follows. We first write: 

.Tec 

where $ is the center manifold function and xj = xj(t) g M are the amplitudes of the critical 
modes. Multiplying the governing evolution equations by the conjugate eigenvectors ipj, J G C 
we get: 

d^J _ « .on™ , (G(^,^),^}) 



(5.7) --f = /3jiR)x 



When the linear part of (5.7) is diagonal, we have the following approximation of the center 
manifold; (see [7]): 

(5.8) -CR^ix,R)=P2Gi^p^)+oi2), 

where Cr = Lji \e2, P2 is the projection onto the span of the critical eigenvectors and 

o(n) = o(N") + 0(N"|/^(i?)|). 

As we have said, we will utilize an expansion of the center manifold using eigenfunctions ej of 
the Laplacian operator: 

s 
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The crucial point here is that P2ej — if J £ C . This is because the eigenvectors of the 
original linear operator and ej's have the same form with different coefficients due to separation 
of variables. So 

(0j,eK>=O, for J ^ if. 

Now, to calculate this approximation we need to compute terms of the form (e^, C-n* es). If 
we denote S = (si, S2, S3), then for S ^C, there are three cases to consider: 

• U si + si = and S3 7^ 0, then 

{esX*Res) = -S37r^(es,es). 



(5.9) 



Using (5.8), the eg coefficient of the center manifold is: 

S^TT^ (65,65) 



• li sl + si 7^ and S3 = 0, then 



(5.10) 



(5.11) 



Using (5.8), the 65 coefficient of the center manifold is: 

*s = — — 2Y + 
P2a|(e|,e|) 

If sf + S2 7^ and S3 7^ 0, then we define: 



S ' 



m,n — 1, 2, 3, 



As can be found by plugging the eigenvectors (5.371 into (5.3): 



Here 



^5 = 



-7i 



P2 7|s37ras^ 



-P2 7s"s^ 



L1L2 



Jo 



o Si7ra;i o S27ra;2 n 

cos cos cos S37ra;3. 

Li L2 



Clearly vs 
nonzero. In this case, denoting 

$5 = 



where n is the number of sub-indices of = (si,S2,S3) which are 



$1 $2 $3 



where $^ denotes the coefficient of 6g, we have: 



(5.12) 



vsAs<^s 



{Gm.e\) (G(r),e|) {Gm.e%) 



By (3.4), we see that for 5 ^ C, 



detA 



Plp27s / D 2 , 2 2^ 2 I 6\ ^ n 

5 |i?=_R,-= ^(-iT^^a^ + S37r Q75 + 75) < 0. 



This ensures that (5.12) can always be solved when S ^ C. 
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5.2. Proof of Theorem [TJ Now let there be a single critical index J = {j, 0, 1) in C. A simple 
calculation yields that 

Thus 
(5.13) 



^{x)=x^ J2 $565 + 0(2), 



where 



Srou = {{0,0,2),{2j,0,0)}. 
To compute the coefficients $5, we multiply ( |5.8[ ) by 65 to get 

{es,C*es) 



(5.14) 



S £Sr 



Thus we can write (5.7) as: 

(5.15) 

where 

(5.16) 



dx 



P{R)x + bx^ + o{3), 



seSr- 



Using (5.9) and (5.101, we find that 
(5.17) 



T_ 
8tt' 



$0,0,2 - 

Also directly computing, we have, 

Gs{tpj,e2j,o,Q,'4'j) 



2j,0,0 



71^72 



4p2a 



2 ■ 



(5.18) 



4a2 



Pl7r272(^2_^2)g^ 



(5.19) 



i^J, rj) = ^^^^(P27'(l + Pi) + Pi(P2 - l)n'Q). 



When p2 > 1, or if p2 < 1 and Q < Qo, (see ( [SJ])) w e see that { t/jj, ip *j) > 0. Thus (^/jj, V'j) 
has no effect on the transition. Next plugging (5.17|-( |5T8 | into (5.16), we find that at the 
critical Rayleigh number, 

(5.20) b . 3,^;^\^y^]^./ 2.^Q(^- - - pW.Rr). 

The transition in this case depends on 6. If 6 < then the transition is Type-I and the 
roll pattern with index J is stable after the transition. On the other hand if 6 > then the 
transition is Type-II and the roll pattern with index J is unstable after the transition. Note 
that b in (5.20) and b in (4.1 ) have the same sign. That finishes the proof of the first assertion. 

Now we will prove the second assertion. For this let J — {j, k, 1) with j ^ 0, k ^ Q and 
consider 

Srec = {(0, 0, 2), (2j, 2k, 0), (2j, 0, 0), (0, 2k, 0), (0, 2k, 2), {2j, 0, 2)}. 
It is easy to check that the center manifold function is: 



(5.21) 



$ = ^ $ses + 0(3). 



seSr 
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As in the proof of the first assertion, the crucial parameter to find is the coefficient of the cubic 
term ( 

(5.22) 



term of the reduced equation which is given in a similar fashion as b in (5.f6) 



By (5.9| and (5.10), we can compute: 



(5.23) 



^"0,0.2 ~ — 



5=(2j,2A:,0),(2j,0,0),(0,2A:,0), 



7 

16n' 



By (5.12) we can compute for S = (2j, 0, 2), (0, 2fc, 2), 



(5.24) 



$1 
$1 



7r72(4Q;2 -al)(R + ns) 
16a27|(i? - i?s) 



7| 

pip2a^iRs + Vs){R + Vs)'^ 
27r 



where 77 and Rs are defined in (4.2 1. To compute a, we also need the following: 

L1L2 



(5.25) 



G's(^,/,e2j,2fc,o>j) 
Gs{ipj,es,tp*j) 



16a 



2Pi7r2Q72(7r2-a2) 



|^Pi^'g7'((-«' + Ic^lW - 5 = (2j, 0, 0), (0, 2k, 0), 



Gs(V'J, 60,0,2,'/'}) = ^^^^Pip27ri?7^, 



For 5 = (2j, 0, 2), (0, 2A:, 2), we have: 
(5.26) 



Gs(V'j,e^,'0}) 



ii-^27i'7^(4a2 - a|) 



P27'* - PlTT^Q 



„4_„.^2, 

pip2i?a2 
-2pi^Q7' 



Now using (5.231, (5.24) and after some simplification, wc find that a in (5.22) has the same 
sign as a in (4.1 1. The rest of the proof is same as in the first assertion. That proves the second 
assertion. 

Finally, notice that by using (3.2), the first critical eigenvalue can be written as /3(-R) = 
c{R — Rc) + o{R — Rc) where c ^ is a constant. That proves the third assertion and the proof 
of the Theorem [1] is finished. 

5.3. Proof of Corollary [1], For the proof we will need the following lemma. 



Lemma 1. Let aj^ be the critical wave number minimizing (3.4). Then 



21/3(22/3 + 1)1/2' 



for all Li, L2,Q. 

Proof. Fix Li,L2> 0. Let 



Ri{x,Q) 



(7r2+x2) 



Due to convexity of i?i with respect to x, for fixed Q, Ri has a unique global minimizer Xj. on 
X > 0. Now, there are two numbers ai, a2 of the form 

'pTT^/Ll + k^n^/Ll, j,keZ 
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such that aj^, — ai or aj^ = a-i where a\ < < a2 and there are no other numbers of this 
form between ai , a2 . 

Since Xr is an increasing function of Q, aj^ is also an increasing function of Q. So it suffices 
to assume that Q = 0. In this case Xr ~ 2.22. We write: 



Jr r2 + r2 ■ 

^1 ^2 



Now fix Li. First, we wih consider that L2/L1 is sufficiently small so that kr = 0. Define 

L{m) = ((m + l)m)i/3((TO + 1)^/3 + m^^^y/'^, to G Z, m > 0. 

We claim that if L{m — 1) < L < L{m) for to > 1 then jr ~ to, i.e. aj^ = In fact, jr G Z 
is an increasing function of L (when fc = 0) and it can be easily checked that 



^1 ', - Ri 



1)tt 



^L{m + l)J \L{m + l) 

Thus, if L(to - 1) < i < L{m) then 

TOTT TOVr TT TT 

" "T - L{mj - i(T) - 21/3(22/3 + 1)1/2 ^ 

Now if L2/L1 is not small, then kr is not necessarily zero and aj^ will be closer to Xr- Thus 
the proof of the lemma is finished. □ 

Now, to prove the Corollary[Tj consider b defined by (4.1 ). We notice that the condition 6 > 
dissects into three conditions: 

(5.27) P2< — , Qn > 2,2(^2 _,2)_p2^2^2- 



To make use of (5.271, it is important to get bounds on the wave number aj. One can see that 
there is no upper bound on the wave number by letting ^1,^2 — ^ 0. However, Lemma [l] shows 
that the minimum wave number can not go below « 1.55. Moreover by definition of Rr, see 
( |3.4[ ), a J will increase as Q increases. And using a similar argument as in the proof of Lemma 
[Ij we find that aj > tt for Q > 306. Thus the transition is Type-I regardless of the container 
size when Q > 306. Next, using the fact that aj > 1.55 and the second condition in (5.271, one 
sees that the transition is always Type-I if p2 > 2.24 regardless of the wave number. 

When at least one of the length scales ^1,^2 is much larger than 1, then the wave number 
will be close to the absolute minimizer of (3.4 1. The relation between Q and the critical wave 
number is then found to be: 

(5.28) Qtt^ = -TT^ + 3n^aj + 2a|}. 

The minimum aj can be found to be aj — tt/\/2 by setting Q ^ above. Then the second 
condition in (5.27) implies that the transition is Type-I if p2 > 2/a/3 > 1.15. Finally by (5.281 
one finds that for Q > 47r^, aj > tt. 

5.4. Proof of Theorem [2J Let / = (7, k, 1) and J ~ (0, 2fc, 1) be the critical indices with 
j,fc>l. 

Consider = (0,0,2), S2 = (2j,2fc,0), S3 = (2j,0,0), ^4 = (0,2fc,0), S5 = (0,4fc,0), 
5-6 = (j,fc,2), ^7 = (j,3fc,2), ^8 = (j,3fc,0), ^9 = (j,/c,0), Sio = (0,2fc,2), Sn = (2j,0,2). In 
this case the center manifold function is: 

^ = a;?(*Siesi + '^S.es,) + xiXj{"Y<S>s^esJ+x^j{<i>-I.^es,+^s,es,). 

4,10,11 7,8,9 
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After letting x = xj, y — xj, the reduced equations (5.7| become: 

dx 

(5.30) 



dt 
dy 



[3x + x[aix'^ + 022/^) + o(3), 



dt 



Py + y{bix' + b2y')+o(?,). 



The coefficients in (5.30) can be found as: 



(5.31) 



ai 



02 



bi 



[$^^G,(V;/,es,,^;)+ 5] $5,G,(7^,,es,,V^;) I , 

y i=2, 3,4,1041 ) 

I $i^G,(V/,es,,^|)+ ^ $s,G.(V^j,es.,^?) I , 

^'"^^^ V i=6,7,8,9 y 

^ I $^s^G.(V'j,es,,V})+ E *s.G.(V7,e5.,V'}) I , 

'"^-^^ V 4=6,7,8,9 / 



(5.32) 



We can find the following relations: 

Gs(-0j,esi,i/'j) = 2G,5(V'/,esi,'(/'|), Gs(i/'j, 65^, -(/;}) = 4G5(V'/, esa, V'/), 
Gs(-07,V'J,ese) = 2G(V'/,V'/,esio), Gs(-0/, V'J, eg^) = 2G(-0/, V/, V'Sn), 
Gs(-07,es8,V'}) = Gs(i/'/,esio,V'/), Gs(V'/, eg^, V'}) = Gs(V'/,esii,'/'7), 

Gs('0J,es,,'0/) = Gs(V'/,es., V'}), « = 6,7,8,9. 



(5.33) 



4 ^ $5^G3(V/,es,,V;) = E $s.G,(V'j,es.,V7*), 



i=3,4 



4 ^ $s.G,(^7,es.,V;)= E <I's.G,(V'j,es.,Vl)- 



2=6,7 



1=10,11 



Using these relations, we can write (5.29) as 

dx 

(5.34) f 



dt 



^ I3x + x {ax^ + 2(2a - b)y'^) + o(3), 
= /3y + y{{2a - b)x^ + 2by'^) + o(3). 



find that a and b are given by (4.1 ) 



Here a = ai and 26 = 62 in (5.31). After simplification of same positive terms appearing, we 



The equations (5.34) have eight straight line orbits in total which are on the lines: x-axis, 



y-axis and x = ±2y. The eight bifurcated solutions of (5.34) are given by: 



(5.35) 



Rec± = ±v/-/3/2fe(0,l), Roll± = ±^J-l3/a{l,0), 
Hex'^ = (-l)V-/3/2(4a- 6) (2,±1), i = l,2. 
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Moreover the Jacobian matrix of the vector field (5.34) has the foUowing eigenvalues at these 
bifurcated points: 



(5.36) 



■^Rec± 

Hex, 



_ 4P{a-b) 

b ' 
_ 4ff(a-b) 
2a ' 
_ 2/3(a-b) 
— 4a-6 ' 



^2 - 
^RoU± — 


~2f3, 


X2 — 


-2/3, 


■^u 1,2 — 
Hex 


-2/3 



Using ( 5.34 1-( 5.36 1, we see that the transition depends on the sign of a, b, Aa — b and a — b. 
These 4 lines separate the a, b plane into eight regions. In each region, it is possible to find 
which solution is bifurcated on /3 > 0, /3 < and the stability of these solutions. The analysis 
is given from Figure [3] to Figure [Tl] These eight cases exhaust all possible transition scenarios. 
The proof is finished. 

5.5. Proof of Corollary [2[ First notice that 







7o,2fe,2 = 4Qg^, 



47r^ < 4a^ - 
- 47r^ < 4a| J. 



47r2 



Thus 



p2 > 8 > 
pi7r'Q(l 
Pi7r2Q(l 



72j",0,2 



4>4m(i 



7} 



1-lL 

2 

72j,0.2 

1^ 
2 

72j,0,2 



) < P1P2 



7! 



42%^)> 
7} 



(7j + T^Q) 



72j.0,2 



(1 



.7|jA2n 



72j,0,2 

4 / 7j 
P27} < Plp2^^ 
72j\0,2 



(7j 



Noticing that Rn < i?2j,o,2 we see that K2j,o < with k defined in (4.2 1. The same argument 
shows also that Ko,2fc < 0. By Lemma [l] tt^ — 5q!j < 0. That shows that a < where a is 
defined in (4.1 ). Now by Corollary [T] 6 < when p2 > 2.24. We can write a/6 as: 



= 1 



> 1. 



b b 
Now, by Theorem |2j the transition depends on a/b and we see that the system is in Region I2 
when a/6 > 1, a < 0, 6 < 0. That proves the corollary. 

5.6. Proof of Theorem [sj Define /3 (R) = 13]^ (R) = X{R)+ ip [R) with A [Re) = 0, p [Re) = 
p > 0. The corresponding eigenvectors and conjugate eigenvectors are 



The horizontal components of the critical eigenvectors can be captured as in (5.2). The ampli- 
tudes of the vertical velocity, the vertical magnetic field and the temperature of ipJa^ V"}^ Eire 
given by: 

Wj^{p^'P + ^j)ajR, 
(5.37) e,7 ^ LojWhl 

H3 J ~ a^jlnR, 



W*j = (/3 + p27})prV 



(5.38) 



HIj = -lirQa] 
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Here ujj{i3) 



and 
(5.39) 



Then 



We define the eigenvectors 



(5.40) = rji - crz, 



Then by (5.391, 



«,<i>}^) = (V'l,cI>}i)=0, 
B:= {^i,'f*j])-' = {ijl,^t) 



^0. 



Now we write = xipj + yipj + $ where $ is the center manifold function, x, y G M. We can 
write the reduced equations of (2.4), (2.5) as 



(5.41) 



dx 

'di 
dy 

dt 



-px + Xy + B{G (V,V'),$}!) 



We have the following approximation formula of the center manifold function, see [7]. 
z = xtpi+y^\ e El, 



For 



where 



{-CrY + Ap' ) {~Cr) a> (z, R) = [{-CrY + 4p^ j P2G (z, i?) 

- 2p2p2G (z, i?) + 2p2p2G (xV-^^ - i?) 
+ p (-£fl) (G {xi^X + - ^V-l, i?)) 

+ G {yi,^ - x^Pl , x^]^ + y^Pl ,R)+o{2), 

oi2) = o[\\z\f)+o{\\{R)\\\zf). 
El ® E2, El — span{^pj^,^j^}, E2 ~ E^ . As in the real case: 



Here Cr = Lr for H 

(5.42) (g (^},, Vl) , V-}) - if J ^ (0, 0, 2) or J ^ (2^,, 2fc„ 0). 
By ( |5.42 ) and the above formula, we have the following approximation for the center manifold 

(5.43) $ (z, R) = $1^-002 + *2V^2j,2fc<=o + o (2) , 
Here $1 and $2 are given by: 

$1 = Aix^ + Asxy + A32/^ 

To compute Ai, we list the necessary results 

(5.44) 5i-:=(G(^}^,Vl),V'c 



-ReLfjReujjS, 
gli = -ImLpReujf), 



9l2 



-ReLpImujjS, 
922 = -ImLpImujp; 
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4:= (g (^},>i)'^2l2.co), 

g^^ = -^Q-^a'jRReKp, gf^ = -2Q-^a\RImKp, 

2 2 2 n 

921 — 512) ^22 — 



cli = p2ReLpReu}-0, c^2 = p2ReLpIniu^, 
C21 = p2l'mL^Reuj, C22 = p2lfnLpImijjj; 



'^ij •= (G'('^/3-^c\V'2i,2fe,o) , 

C21 = ^ajJmKjs, = 0; 



rfii = 2(a^^ - '!r^)ReK^, di2 = {a^j^ - Tr'^)ImK-^, 
d2i = {aj^ - ir'^)ImK0, ^22 = 0. 



5p = P + p2lX, 

K0 = -^QRaln%%, 

L0 = ^llalRnp^'Sp. 
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Using these resuhs, the coefficients of the center manifold function can be computed as foUows: 



A2 
A3 

A5 

A. 



47r2 (167r4 + 4p2) 

{{gh + gli) IStt^ + (gii - .92^2) (p4^^ + 1)) 
47r2 (167r4 + 4p2) 

(,gi2 (16^4 ^ 2p2) + 5i^2p2 + giip47r2 + .gi^) 



(I6p 



2 4 



47r2 (167r4 + 4p2) 
2p2) _ (4pp2al + 1) 9I2 



4p2aj^ (I6p 



2^,4 



V) 



(4pp2«l + 1) gli + 32p2a5^g2_^ 
4p2a2^(l6p2a4^+4p2) ' 

2p'g?i + (4pp2al + 1) ff?2 
4p2a2^(l6p2a4^+4p2) 



Plugging! 5.43 1 into (5.411, we obtain the following ODE: 

^ = Ax + + 0302;^ + a\^x^y + al^xy^ + al^y^ + o(3), 

^ = ~px + \y + aj^x^ + al^x^y + aj^xy'^ + al^y-^ + o(3). 
Using the approximation ( 5.43[ ), we get 



(5.49) 



(5.50) 



(5.51) 



(G (V, V') , = x'^i - Cc\^) + a;<i>2 - Cc?^) + 

a;$2 (dii - Cdi2) + 2/$2 (^21 - 6*^22) + o(3), 
(G (V, V-) , 'f!/!) = ^i;*! {Cc\^ + c\^) + a;$2 {Ccl, + c?^) + 



(Gc^i + 0^2) ■ 
a;$2 (Gdii + ^12) 



C21 + C22) 



{Ccl 

-y$2 (Gd2i +^22) + 0(3). 



Now we give the coefHcients of the reduced equation (5.49). Using (5.43), (5.50) and (5.51) 



the coefficients are computed by plugging (|5.50|), (|5.5l|) into (|5.4l|) and are given by: 



^30 



^21 



^12 



^21 



^12 



S(Ai(cii-wH2 
BiAiicl, - Ccl^ 
B[A2[c\^^Cc\^ 
BiAsicl, ~ Ccl^ 
BiA,iCcl^ + c\^ 
B{A,{Ccl, 
B{A2{Cc\, 
B{A^{Cc\^ 



22 



22 



X 

Y = c 



12 



il2 



11 



+ A4X), 

+ yl2(c}i - Cc\^) + A^{cl^ + d2l) + A^), 
+ A3(c}i - Gc}2) + A5(C21 + dil) + AgX), 
+ A6(c2i +d2l)), 

+ A2(Gcii + c}2) + A^C{cl^ + rf2i) + A^Y, 
+ A3(Gc}i + c}2) + A^C{cl^ + rf2i) + A^Y, 
+ A6G(c2i +d2l)), 

-Cfl2), 

dii). 
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The transition of (2.4)-(2.5) is determined by the sign of the following number at R 
see [7], 



(5.52) 



Stt 



^30 



2 \ /I 



•■21) 



which has the same sign as b defined by: 

Di + D2 



(5.53) 
Here 



Qtt 



nRr 



7r2 (167r4 + V) 2p27}^ {^^PI<^X 



Di = 2p2 (37,1^1 + + tI^s) (i^iV'ii + £^2^-21) , 
D2 = 2p2 (Ml + 7,1^2 + SpAa) (£;iV2i - £^2^-11) , 

= ^4(3^-11 + 2^21PP2^X,'^) + ^5'/'21 + ^6(V'll + 2V'21/OprS7/). 

E2 



(p2+pi)(7lpr' + Q^'(pi + i)"'), 

(P2pl)"^(p2 +Pl)p7l, 

- ((I67r'^ + 2p2) ^_2^i;i + 2p^V£^2P^ - HiE2pT^'^ - Eip^^p) 

- (16^4 (i?27,l + P2'PE1) + (7,1^1 - P2~ V^2) {W + 1)) : 

- {p2^pE2 (16^4 + 2p2) + 2p2f^^iJ, + 47r2p- V'i?i + 7li?2) 



A2 
A3 

A4 - 27I (I6p^a*,^ + 2p2) _ (4pp2a2^ + l) p. 
A, 
A, 

^^21 = (prVai^c + p2i^i-B2)7; 



27}^ (4pp2a',^ + 1) + 32p2p, 
V7I + (4pp2a^^ + 1) p, 
= -((p2pi)-V'«l^c + p2i?|7c), 



Finally from (3.5), we notice that 
(5.54) lX=0[Q 



)l/3 



as Q 



Using this, we see that, as Q — )■ 00, 



(5.55) 



Pip2(l-p2)7r^ 

pl + 1 



Q. 



for c > 0. Plugging (5.551 into the expression h defined in (5. 53), we see that 5 < as Q — > 00. 
Also, it can be shown that the same result holds in the limit of small p. Theorem [3] is proved. 

6. Physical Remarks and Conclusions 

In this work, we investigate several transition scenarios of the magnetohydrodynamics (MHD) 
equations. As is well known, for (MHD) due to non-selfadjoint linear opeartor, the transition 
can be caused by a finite set of real or nonreal eigenvalues crossing zero. 

When the first eigenvalue is real and simple, the transition depends on the character of the 
first critical index. In this case, the transition can only be Type-I or Type-II depending on 
a number exactly given by (4.1 1. In particular, when the first critical eigenmode has a roll 
structure, the type of transition is independent of the Prandtl number pi. In this case, the 
transition is Type-I if Q > or p2 > p* where and p* depends on the length scales of the 
domain. We find that Q* < 307 and p* < 2.24 regardless of the length scales of the domain. As 
max{Li,L2} t 00, Q* i 47r2 and p* | 2/\/3. 
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Next we study the case where there are two critical real eigenvalues. In this case we only 
consider the special geometry 

Li ^ J 
L2 kVi' 

with positive integers j, k and Li, L2 denoting the horizontal length scales of the box. With 
this assumption, it is possible that two modes which can characterize a hexagon pattern become 
unstable at the same critical parameter. In this case we find that all types of transitions are 
possible in a total of eight different transition scenarios. However, in our numerical investigation, 
we encountered only two of these scenarios. We find that when Q < Q* a-^d p2 < P* the system 
moves from a Type-Ill transition regime to a Type-I regime as p2 crosses p*. The minimal 
attractors in the stable domain of the Type-Ill transition regime have a steady rectangular 
pattern. In the Type-I transition regime, the rectangles and rolls are minimal attractors, and 
hexagons are unstable patterns after the transition. In this case we prove that for p2 > 8, the 
transition is always Type-I with rolls and rectangles as stable patterns and hexagons as unstable 
patterns after the transition. However p2 > 8 is a crude estimate and our numerical investigation 
suggests that this type of Type-I transition will be preferred for p2 > p2 where p2 < 2.24. 

Finally, we consider the case where the first eigenvalue is simple and nonreal. This is always 
the case when p2 < 1 and Q > Qo- We only consider a roll type critical eigenmode. In this 
case the transition can be Type-I or Type-II. In particular for Q sufficiently large or when the 
oscillation frequency p is sufficiently small, the transition is Type-I and the transition structure 
is a time periodic roll pattern. 
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